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Abstract 

We present a model for the description of nuclear matter and finite nuclei, and at the same 
time, for the study of medium modifications of nucleon properties. The nucleons are described 
as nontopological solitons which interact through the self-consistent exchange of scalar and vector 
mesons. The model explicitly incorporates quark degrees of freedom into nuclear many-body 
systems and provides satisfactory results on the nuclear properties. The present model predicts a 
significant increase of the nucleon radius at normal nuclear matter density. It is very interesting 
to see the nucleon properties change from the nuclear surface to the nuclear interior. 
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I. INTRODUCTION 



One of the most interesting topics in nuclear physics is to study how the nucleon prop- 
erties change in nuclear medium. So far, there are many experimental evidences indicating 
that the properties of the nucleon bound in nuclei are significantly modified from those of 
a free nucleon. The famous European Muon Collaboration (EMC) effect shows that the 
nucleon structure functions in nuclei deviate from those in a free nucleon [l|. Important 
evidences for medium modifications also come from recent polarization transfer experiments 
at the Thomas Jefferson National Accelerator Facility, which observed a difference in the 
electromagnetic form factors of a proton bound in a helium nucleus compared to a free 
one Q . On the other hand, there are numerous theoretical works on the study of in- medium 
nucleon properties based on various models 0, \i, ^,0]. At present, we are still far away from 
describing nucleons and nuclei in terms of quarks and gluons using quantum chromodynam- 
ics (QCD), which is believed to be the fundamental theory of strong interactions. Hence it 
is highly desirable to build models which could incorporate quark degrees of freedom and 
respect the established theories based on hadronic degrees of freedom. 

The quark-meson coupling (QMC) model proposed by Guichon 0] can be considered 
as an extension of the extremely successful theoretical treatment of nuclear many-body 
systems, known as quantum hadrodynamics [8], to include the internal structure of the 
nucleon. The QMC model describes the nuclear system as nonoverlapping MIT bags in 
which the confined quarks interact through the self-consistent exchange of scalar and vector 
mesons in the mean-field approximation. In the QMC model, the quark structure of the 
nucleon plays a crucial role in the description of nuclear matter and finite nuclei jf|. It is 
also possible to study the medium modification of nucleon properties in the QMC model 
where the quark degrees of freedom are incorporated explicitly in the nuclear many-body 
system jj. In the past decade, the QMC model has been extensively developed and applied 
with reasonable success to various nuclear phenomena 0,0, HQ, HQ HQ. There 



are also other models that incorporate quark degrees of freedom in the study of nuclear 
many-body systems. The quark mean-field (QMF) model 0] takes the constituent quark 
model for the nucleon instead of the MIT bag model, where the constituent quarks interact 
with the meson fields created by other nucleons. The QMF model has been successfully used 
for the description of nuclear matter, finite nuclei, and hypernuclei Il9l |. Recently, 

the QMF model has been extended to a model based on SU(3)l x SU(3)r symmetry and 
scale invariance Q]. Using the Nambu-Jona-Lasinio model to describe the nucleon as a 
quark-diquark state, it is also possible to discuss the stability of nuclear matter based on the 
QMC idea 21[. The main advantage of these models is their simplicity and self-consistency 
in taking into account quark degrees of freedom in the study of nuclear many-body systems. 

In this paper, we take the nontopological soliton bag model originally proposed by Fried- 
berg and Lee 22J], which is also called Friedberg-Lee model in the literature, for the descrip- 
tion of nucleons in nuclear medium. In the Friedberg-Lee model, the nucleon is described 
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as a bound state of three quarks in a nontopological soliton formed by a scalar field with 
nonlinear self-interactions. The Friedberg-Lee model has the benefits that it is manifestly 
covariant and it exhibits dynamical bag formation due to the coupling of quarks to the 
phenomenological scalar field. Furthermore, it includes the MIT bag as a special case. The 
soliton bag model has been extensively used to study the properties and structure of hadrons 
in free space 23, 24, 2s|, 26, 27, 28, 29l[30, 31]. It has also been applied to discuss the medium 
modification of nucleon properties [3j, |32| and to study the dense matter properties within 
the Wigner-Seitz approximation In the present work, we develop a model to study the 
properties of nuclear matter and finite nuclei by describing the nuclear many-body system 
as a collection of nontopological soliton bags. The quarks inside the soliton bag couple not 
only to the scalar field that binds the quarks together into nucleons, but also to additional 
meson fields generated by the nuclear environment. The nucleons interact through the self- 
consistent exchange of these mesons treated as classical fields in the spirit of the QMC and 
QMF models. Because in this model the soliton solution is significantly changed by the 
additional meson fields, it is possible to investigate the modification of nucleon properties 
in nuclear medium. 

This paper is arranged as follows. In Sec. [Til we briefly describe the Friedberg-Lee model 
for nucleons both in free space and in nuclear medium, and then discuss the properties 
of nuclear matter as a collection of nontopological soliton bags. In Sec. IIIIl we show the 
properties of finite nuclei with the nucleons described as soliton bags, and also discuss how 
the nucleon properties change inside nuclei. Section [TV] is devoted to a summary. 



II. NUCLEAR MATTER AS NONTOPOLOGICAL SOLITON BAGS 

In this section, we first give a brief description of the Friedberg-Lee model for individual 
nucleons. Then we develop a model for nuclear matter regarded as a collection of nontopo- 
logical soliton bags. This model enables us to investigate possible modifications of nucleon 
properties in nuclear medium. 



A. Free nucleon 

In the Friedberg-Lee model, a single nucleon is described as a bound state of three quarks 
in a nontopological soliton formed by a phenomenological scalar field with nonlinear self- 
interactions. The Friedberg-Lee model in its simplest form is implemented through the 
effective Lagrangian density 

£ = ^{i lfl d» - m - + X -d^d^ - U(<f>), (1) 

where ip denotes the quark field. The quark mass m is usually taken to be zero for u and 
d quarks. is a color-singlet scalar field that may be interpreted as the phenomenological 
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representation of quantum excitations of the self-interacting gluon field. The self-interaction 
of the scalar soliton field is described by the potential 

W) = |0 2 + + ^ + B. (2) 

Here, the polynomial terminates in fourth order to ensure renormalizability. The constants 
a, b, and c are fixed within a range so that U ((f)) has a local minimum at = and a global 
minimum at = <f> v . The constant B is determined to make U(<f> v ) = 0, and then the value 
£7(0) = B is to be identified with the bag constant or volume energy density of a cavity. We 
note that 0„ is the value of the soliton field in the physical vacuum where the quark gets a 
mass ~ g<f) v from its coupling with the soliton field. Inside the nucleon where valence quarks 
exist, the soliton field is reduced to be near zero that means the perturbative vacuum 
is restored. It is energetically favorable for three quarks to be localized in a cavity in the 
soliton field, referred to as the soliton bag, so that the nucleon appears as a bubble in the 
vacuum. In the mean-field approximation, the soliton field is treated as a classical field 
which is a time-independent c-number field 0(r). The quark field operator is expanded in 

a complete orthogonal set of Dirac spinor functions as ip = J2bk?Pk, where bk is the fermion 

k 

annihilation operator. For a nucleon, the three valence quarks are in the lowest Dirac state 
■00, then and ip satisfy the coupled differential equations 

(-id- V + fi'/30)^o = eo^o, (3) 

- V 2 + ^ = -3# ^o. (4) 

The coupled equations have to be solved numerically. For the Dirac wave function of the 
lsi/2 quark, we use the notation 

*> = f .-"i* , , | X, (5) 
\ia-rv (r) 

with x — I 1 or ( ) being the Pauli spinor. The total energy of the nucleon is given by 







E = 3e + 47r J drr 2 \ (^j + U(<t>) 



(6) 



The mean-square charge radius of the proton is given by 

(r 2 ) =4n J dr r 4 (u 2 + v 2 ) , (7) 

and the proton magnetic moment is given by 

87T r , , /o , 

li p = — J dr r uv. (8) 

The ratio of the axial-vector to vector coupling constants is given by 

9a/9v = ^Jdrr 2 (u 2 - ^ . (9) 
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B. Center-of-mass correction 



There are several methods for calcula ting the center-of-mass (cm.) correction to nucleon 



properties in the soliton bag model 24], |25|, |26j. In this paper, we consider two approaches 
for the cm. correction so that we can estimate how sensitive the results are to the cm. 
correction approach used. First, we adopt an operationally simpler approach based on the 
relativistic energy-momentum relation to take into account the cm. correction, which has 
been extensively discussed in Ref. [241] • The rest mass of the nucleon in this approach is 
given by 

'E 2 - (P 2 ), (10) 



M 



where P is the total-momentum operator. The corrected root-mean-squared (rms) radius is 
given by 



\ 



2eo 3e| 
E E 2 



(r 2 ) + 



2E r 



(11) 



We also use the Peierls-Yoccoz projection technique to evaluate the cm. correction to 
nucleon properties. Following the method described in Ref. 25j, the rest mass of the nucleon 



is given by the expectation value of the Hamiltonian in the zero-momentum projected state. 
In this approach, the nucleon state is assumed to be the direct product of a soliton coherent 
state and a three-quark state 

\N)=e XAi b\blbl\0). (12) 
The soliton coherent state is given by 



e XA '\Q) =TJexp 



/k4 



|o>, 



(13) 



where is the creation operator of the soliton field. The soliton field operator can be 
expanded in terms of and 



+ (2vr) 



-3/2 



/ ^1 f 



a k e ikr + 4e~ ikr 



(14) 



The expectation value of the soliton field in the coherent state has exactly the properties 
of the mean- field value <p obtained by solving Eqs. ([3]) and (T4j). The soliton coherent state 
is localized, so the nucleon state has no definite momentum. To construct a momentum 
eigenstate, we use the Peierls-Yoccoz projection. The zero-momentum projected state is 
given by 

|P = 0) = J d 3 X\X.}, (15) 
where |X) is a nucleon state localized at the point X 



|X)=e^ t ( x )6 t 1 (X)4(X)4(X)|0). 



(16) 



5 



The operator A^X) can be written as 

AAt(X) = / d'k^f^al 



(17) 



For a translationally invariant operator O, its expectation value in the zero- momentum 
projected state is given by 



(O) 



(P = 0|O|P = 0) _ Jd 3 Xd 3 Y(X\Q\Y) _ $ d 3 Z(-\Z\Q\\Z) 
(P = 0|P = 0) ~ Jd 3 Xd 3 Y(X\Y) ~ jd 3 Z{-\Z\\Z) 



(18) 



where Z = Y — X. The normalization condition can be expressed as 

(_I Z |iz)=iV (Z)iV g (Z) 3 , 



where 



N^Z) = exp 



N q (Z) = J d 3 r^ (r+iz) ^ (r-^ 



1, 
2' 

Z). 



(19) 



(20) 
(21) 



Therefore, the nucleon mass is given by the expectation value of the Hamiltonian in the 
zero-momentum projected state 

M={:H:) = {:H q + H q ^ + H <t> :). (22) 

These expectation values of normal-ordered products are 

Jd 3 r4 (r + |Z) [-ia ■ V + g/ty(r; Z)]jj (r - |Z 



(: H q + H q(j) :) = 3/ d^ZN^N^ 2 

= 



Jd 3 ZiV (Z)iV g (Z) 3 
!d^ZN 4> {Z)N q {Zf£ 4> {Z) 



J#ZiV^(Z)iV g (Z) 3 



where 



2^( r ; z ) 



i|V0(r;Z)| 2 + f/(0) 



The expectation values of the soliton field operators are given by 



(r;Z) 



r-iz 

2 



7r(r;Z) = -z (2tt) 



-3/2 



2 



/ k ( -Z ) e^ r - / k 



-~Z ) e" ik - r 



(23) 
(24) 

(25) 

(26) 
(27) 



We note that the corrected rms radius in this approach has the same expression as given by 
Eq. flnj H- 
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C. Parameters in the Friedberg-Lee model 

In the Friedberg-Lee model, the parameters a, b, c, and g are constrained by reproducing 
reasonable nucleon properties. We adjust the parameters to fit the nucleon mass M = 939 
MeV and rms radius r c = 0.83 fm. In the present work, we take two sets of parameters using 
the cm. correction given by Eq. (fTOj) . Set A: a = 0, b = —79.61 fm -1 , c = 780, g = 13.7 is 
characterized by a = 0, where U(4>) has a inflection point at = 0. Set B: a = 69.945 fm~ 2 , 
b = —1600 fm -1 , c = 12200, g = 24.55 is characterized by B = 0, where the relationship 
among the parameters b 2 = 3ac is obtained by requiring B = 0. The parameters in the 
Friedberg-Lee model have to be in a range where U(4>) has a local minimum at <p — and 
a global minimum at = <f> v . Therefore, sets A and B correspond to the two limiting cases 
that can be seen in Fig. [TJ 

In order to compare the results with different cm. correction approaches, we take another 
parameter set using the projection technique for the cm. correction given by Eq. (1221) . Set 
C: a = 19.08 fm" 2 , b = -335.02 fm" 1 , c = 1961.0, g = 24.55 is also characterized by B = 0. 
We use the same coupling constant g in sets B and C. By comparing the results with these 
two parameter sets, we can estimate how sensitive the results are to the cm. correction 
approach used. 

We compute nucleon properties using these three parameter sets. The nucleon mass 
M = 939 MeV and rms radius r c = 0.83 fm are obtained because they are constrains for 
the parameter sets. Set A gives the proton magnetic moment fi p = 2.80 and the ratio 
of the axial- vector to vector coupling constants gA/gv — 0.87. Set B predicts \i p = 2.77 
and gA/gv = 0.90, whereas set C gives \i p = 2.85 and gA/gv = 0.80. We note that the 
experimental values are \i v = 2.79 and gA/gv — 1-25. It is shown that these three parameter 
sets in the Friedberg-Lee model can give reasonable results for nucleon properties in free 
space. 

D. Nuclear matter 

We develop a model for nuclear many-body system based on the Friedberg-Lee model in 
the spirit of the QMC model. Nuclear matter in this model is considered as a collection 
of nontopological soliton bags. The solitons interact through the self-consistent exchange 
of er, u, and p mesons that are treated as classical fields in the mean-field approximation. 
For a soliton embedded in nuclear matter, the quarks inside the nucleon couple not only to 
the soliton field <p which binds the quarks together into the nucleon, but also to additional 
meson fields a, u, and p generated by other nucleons in nuclear medium. We assume that 
the meson mean fields a, u, and p can be regarded as constants in uniform matter and the 
soliton field that serves to bind the quarks together does not participate in nucleon-nucleon 
interactions. Therefore, <fi depends on spatial coordinates inside a nucleon, whereas a, u>, 
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and p are constants. With the presence of these additional meson fields in nuclear matter, 
the quark and soliton fields in the nucleon satisfy the coupled equations 



(-id ■ V + g(3(j) + g q a (3a + g*u + g q p r 3 pj tp = eo^o, 



- V 2 6 + 



dU{4>) 



-3gifj ipo, 



(28) 
(29) 



where g%, g q , and g q p are the coupling constants of the a, u, and p mesons with quarks, 
respectively. We solve the coupled equations and calculate the in-medium nucleon properties 
analogously to the case of free nucleons. Here, the constant a field provides an additional 
scalar potential to the quarks and as a consequence changes the solutions ipo and <fi from 
those obtained by Eqs. (j3J) and (J4]). On the other hand, the uj and p fields do not cause 
any changes of ipo an d 4> except to shift the energy level by a constant vector potential, 
e (a, u, p) = eo (cr) + g^ui + g q p r 3 p. Hence, ipo and can be expressed as a function of the a 
mean field. 

Analogously to the case of free nucleons, we use two approaches to take into account the 
cm. correction for nucleons in nuclear matter. The first approach is based on the relativistic 
energy-momentum relation 33J, in which the effective nucleon mass is given by 



M* 



o- 



E(a) 



where 



E {a) = 3e (a) + An J dr r 2 




+ U{<f>) 



(30) 



(31) 



The second approach is based on the Peierls-Yoccoz projection technique as described in 
Sec. Ill Bt The effective nucleon mass is given by the expectation value of the Hamiltonian 
in the zero-momentum projected state 



M* (a) = (:H:) 



(P = Q| : H : |P = 0) 
(P = 0|P = 0) 



(32) 



which can be calculated analogously to the case of free nucleons. We note that the quark 
wave function ipo and the soliton field <fi are altered by the a mean field in nuclear matter. 
Therefore, the calculated nucleon properties are different from those in free space. They can 
be expressed as functions of the a mean field. The density dependence of these quantities 
is obtained by a self-consistent determination of cr at a given nuclear matter density. The 
}. ( 1301 has also been used in the QMC and QMF models for removing the 
13l . Il7| . Another prescription used in the QMC model is to incorporate 
the cm. correction into the parameter z and assume that z is independent of the matter 
density jllj ]. 

It is interesting to compare the two cm. correction approaches used in this paper and 
investigate how the cm. correction changes in nuclear matter. We define the cm. energy 



prescription o 
cm. motion 
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of a nucleoli in nuclear matter as E c , m . = E — M*, where M* is given by Eq. (1301 in the 
first cm. correction approach and by Eq. f[3"2"j) in the second one. In Fig. [21 we show E cm ^ 
as a function of the effective quark mass m* = g^a, which is proportional to the a mean 
field. The results of the first cm. correction approach are shown by the solid lines, while 
those of the second approach are shown by the dashed lines. We have m* = in free space 
and m* ~ —240 MeV at normal nuclear matter density. It is seen that the variation of 
E cm ^ in nuclear matter depends on the cm. correction approach and the parameter set. 
As the density increases, -E c .m. obtained in the first cm. correction approach (solid lines) 
decreases in the upper panel and increases in the lower panel, whereas the results of the 
second approach (dashed lines) slowly decreases in both panels of Fig. [2j We find that the 
tendency of i? c . m . with the second cm. correction approach is quite similar to that shown in 
Fig. A.l of Ref. [llj which is the results of a serious calculation for a relativistic harmonic 
oscillator potential. Therefore, the treatment for the cm. correction by the projection 
technique is considered to be more reliable than the simple cm. correction based on the 
relativistic energy-momentum relation. 

To investigate in-medium nucleon properties and nuclear matter characteristics, we take a 
hybrid treatment for nuclear matter. The effective nucleon mass and couplings are obtained 
at the quark level, whereas the nucleon Fermi motion is treated at the hadron level. To 
perform a many-body calculation for nuclear matter, we start from the effective Lagrangian 
at the hadron level within the mean-field approximation 



-,mla 2 + -mlu 2 + -m 2 p p 2 , (33) 



£rmf = ^ - M* (cr) - g^u - g p Tr 3 p\ i/> 

1 2 2,1 22,1 

— m„a H — m, .00 H — 

2 2 w 2 

where tp denotes the nucleon field. The effective nucleon mass M* (cr) is obtained in the 
Friedberg-Lee model, which has been given by Eq. ( 1301) in the first cm. correction approach 
and by Eq. (1321) in the second one. The nonvanishing meson fields are replaced by their 
expectation values cr = (cr), u = p = (p 03 ) that are constants in a static infinite 

nuclear matter. The nucleon-meson couplings are related to the quark-meson couplings as 
guj = 3g2, an d g P = g q p [171 ] . From the Lagrangian given by Eq. ( |33l) . we obtain the equations 
of motion for nucleons and mesons in nuclear matter, 

[i7 M 3" - M* (cr) - g u >fu - g P i°r z p] if, = 0, (34) 

2 dM*(a) /Tl . 
m 2 a a = ^-^(^), (35) 

mlu = g w (^°ip), (36) 

m lp = g P (*Pi 0r 3^)- (37) 

With M* (cr) obtained at the quark level, we solve the coupled equations self-consistently, 
and then calculate the nuclear matter properties and medium modifications of nucleon prop- 
erties. In the present model, the quark-meson couplings g%., g%, and g q are determined by 
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reproducing the nuclear matter equilibrium density (0.15 fm~ 3 ), energy per nucleon (—16 
MeV), and symmetry energy (35 MeV). The meson masses are taken to be m a = 500 MeV, 
m oj — 783 MeV, and m p = 770 MeV. We list in Table [I] the resulting nuclear matter prop- 
erties corresponding to the parameter sets A, B, and C used in the Friedberg-Lee model. 
It is shown that the present model can provide a satisfactory description of nuclear matter 
properties. 

Having determined the a mean field self-consistently in nuclear matter, we investigate 
in-medium nucleon properties with the quark wave function obtained by solving Eqs. (|28|) 
and (129]) . In Fig. [3] we plot the Dirac solutions {u, v) and the soliton field (0) as functions of 
the nucleon radius (r) for the three parameter sets used. The results of the nucleon in free 
space are denoted by uo, vo, 4>o, while those at the density p an d 2p are denoted by u\, v\, 
(pi and U2, V2, <p2, respectively. Here po = 0.15 fm -3 is the normal nuclear matter density. 
By comparing the results in medium with those in free space, we find that the quark wave 
functions in nuclear matter are significantly different from those in free space. This leads to 
the medium modifications of nucleon properties because they are calculated using the quark 
wave functions. In the present model, the quark wave function depends on the a mean 
field, hence the density dependence of nucleon properties is obtained through the density 
dependence of the a mean field. It is interesting to compare the results obtained with the 
three parameter sets used. Set A gives stronger in-medium modifications than sets B and 
C, especially at larger r. The results of set B is quite similar to those of set C because both 
of them are characterized by B = 0. It is seen in Fig. [3] that the soliton field ((f)) around 
the surface at high density in the case of set A is much broader than that in set B. This 
yields more leakage of the quark wave function in set A than in set B. It may lead to the 
difference in the density dependence of nucleon properties, such as r c and Qa/qvi because 
they are sensitive to the behavior of the quark wave function at large r. By comparing with 
the sharp boundary of the MIT bag model, we find that set B is closer to the MIT bag model 
than set A. Therefore, the results of set B are rather close to those of the MIT bag model. 
This indicates that the surface behavior of the confinement potential plays an important 
role in determining nucleon properties. 

In Fig. HI we present the ratio of the effective nucleon mass in nuclear matter to that in 
free space, M*/M, as a function of nuclear matter density p. The solid and dashed lines 
correspond to the cases of set A and set B where the effective nucleon mass M* is given 
by Eq. ( 1301) using the first cm. correction approach. The dotted line corresponds to the 
case of set C where M* is given by Eq. ( 1321) using the second cm. correction approach. 
We note that the effective nucleon mass in the present model is calculated at the quark 
level, which is not a simple linear function of a as given in the Walecka model s|. It is 
more like the characteristics of the QMC model. As shown in Fig. HI the effective nucleon 
mass decreases with increasing density, and the results depend on the parameter set used. 
The large difference between sets B and C is mainly due to the different cm. correction 
approaches used in these two cases. The first cm. correction approach used in the case of 
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set B provides a strong increase of E c , m . in nuclear matter (see solid line in the lower panel 
of Fig. [2]), whereas the second cm. correction approach used in the case of set C gives a 
weak decrease of E cm (see dashed line in Fig. [2]). This leads to the difference in M* between 
sets B and C due to the relation M* = E — E c _ m _. Therefore, the drop of M* in set B is 
much more rapid than that in set C. On the other hand, the difference between sets A and 
B can be understood as a consequence of the different behaviors of E c _ m _ in these two cases 
(compare solid lines in the upper and lower panels of Fig. [2]), though the first cm. correction 
approach is used in both sets A and B. 

We show in Fig. [5] the ratio of the nucleon rms radius in nuclear matter to that in free 
space, r*/r c , as a function of nuclear matter density p. It is very interesting to see the 
expansion of the nucleon size in medium. We find the nucleon rms radius increases by 
about 10-16% at normal nuclear matter density. This result is quite different from those 
obtained in other models. For example, the QMC model predicts only 1-3% enhancement 



in the nucleon rms radius at normal nuclear matter density [10], and the QMF model gives 
about 5-9% increase [17|. The chiral quark-soliton model predicts a 2.4% enhancement, while 
the swelling constrained by quasielastic inclusive electron-nucleus scattering data is less than 
6% 34j. We note that although the Friedberg-Lee model can give similar quark distributions 
and nucleon properties to those obtained in the MIT bag model in free space, the quarks 
satisfy different equations and boundary conditions in these two models, and therefore the 
swelling of the nucleon rms radius could be quite different between the QMC model and 
the present calculation. In Fig. [61 we present the ratio of the proton magnetic moment in 
nuclear matter to that in free space, //*/ p p , as a function of nuclear matter density p. It is 
shown that the results depend on the parameter set used. The difference between different 
parameter sets increases with increasing density. Set A gives nearly the same low-density 
behavior as set B, but rather different results at high density. The enhancement in set C is 
smaller than those in sets A and B. We show in Fig. [7] the ratio of the axial- vector to vector 
coupling constants in nuclear matter (9a/ 9v)* to that in free space 9 a/ 9v as a function 
of nuclear matter density p. The results of set A decrease slightly at lower densities, and 
then increase at higher densities. On the other hand, the results of set B drop significantly 
with increasing density. The results of set C is between sets A and B. At normal nuclear 
matter density, we obtain (g A / ' g v )* / \g A / ' g v ) ~ 0.98 for set A, {g A /g v )* l{gAlgv) - 0.87 
for set B, and (9 a/ 9vY I(9a/ 9v) — 0.94 for set C. It is obvious that medium modifications 
of nucleon properties in this work depend on the parameter set used in the Friedberg-Lee 
model. Because set A gives stronger in-medium modifications of the quark wave functions 
than sets B and C as shown in Fig. [3], we obtain larger increases in r* and p* with set A, as 
shown in Figs. [5] and [61 On the other hand, the difference in (9a/ 9v)* shown in Fig. [7| could 
be due to the competition between changes of u and v according to Eq. ([9]). 
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III. PROPERTIES OF FINITE NUCLEI 



In this section, we extend the present model to study the properties of finite nuclei and 
the modification of nucleon properties in a nucleus. The nucleus is described as a collection 
of nontopological soliton bags that interact through the self-consistent exchange of a, ui, and 
p mesons. In principle, these meson mean fields are functions of the spatial coordinates in 
the nucleus, but it is rather complicated if the variation of these quantities over the small 
nucleon volume is taken into account. Therefore, we take some suitably averaged form for 
the meson mean fields in order to make the numerical solution feasible. We use the local 
density approximation which replace the meson mean fields by their value at the center 
of the nucleon and neglect the spatial variation of the mean fields over the small nucleon 
volume U Q 0- The equations of motion for nudeons and mesons in a spherically 
symmetric nucleus are given by 



irtpP - M* (a) - g^cu - g pl °r 3 p - e^-^-j°A 



V> = 0, (38) 



(-A + mj)*=-®^<W>, (39) 
(-A + ml)tj = g u (^), (40) 
(-A + m 2 p )p = g p (4, 1 T 3 iP), (41) 

-AA = e <^pVV>, (42) 

where the mean fields are functions of the radial coordinate of the nucleon center in the 
nucleus. We solve the preceding equations self-consistently with the effective nucleon mass 
obtained at the quark level. 

We present the numerical results of several spherical nuclei. In Table HIl the calculated 
bindin g en ergies per nucleon and rms charge radii are compared with the experimental 



values 35j. By comparing the results of sets A, B, and C, we find that the binding energies 
of set C are larger than those of sets A and B. This is mainly because set C gives smaller 
incompressibility (K = 184 MeV) than sets A (K = 302 MeV) and B (K = 384 MeV) 
when they have the same saturation density (po — 0.15 fm -3 ) and energy per particle 
(E/A = —16 MeV) as shown in Table [B We list the calculated single-particle energies 
and spin-orbit splittings in Tables II I II and \TV\ and compare with the experimental data 
taken from Ref. [jifjj]. It is seen that all calculated spin-orbit splittings are smaller than 
the experimental data, while set C gives the smallest values among the three parameter 
sets. This is due to the large effective nucleon masses obtained in the present model as 
shown in Table [B According to analysis with many quantum-hadrodynamics models, there 
exists a tight correlation between the spin-orbit splitting of finite nuclei and the effective 
nucleon mass in nuclear matter at saturation density, and the spin-orbit splitting increases 
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with decreasing M*. It is known that M*/M ~ 0.6 is required in order to reproduce the 
empirical spin-orbit splittings of finite nuclei. This is the reason why we get small spin-orbit 
splittings in the present model and the smallest values in set C. This shortage might be 
improved by including nonlinear meson self-interactions as taken in the QMC model JjJ. 
We plot in Fig. [8] the resulting charge density distributions for 208 Pb and compare with the 
experimental values |37|. As seen in this figure, the calculated results are in good agreement 
with the experimental values, and there is no explicit difference between the three parameter 
sets. 

It is also possible to investigate the modification of nucleon properties in finite nuclei. 
Using the local density approximation, the nucleon properties at the radial coordinate r in 
a nucleus are obtained through the values of cr(r), because the nucleon properties in this 
model are functions of the a mean field obtained at the quark level. In Figs. [9] and [TDJ we 
show the ratios of the proton rms radius and magnetic moment in 208 Pb to those in free 
space as functions of the radius r. It is found that the proton radius and magnetic moment 
increase significantly at the center of 208 Pb. These quantities decrease to the values in free 
space from the center to the surface of the nucleus. It is shown that the results depend 
on the parameter set used, which is consistent with the results of nuclear matter shown in 
Figs. [5] and El 



IV. CONCLUSION 

By treating the nucleons as nontopological soliton bags, we have proposed a model for the 
description of nuclear matter and finite nuclei, and at the same time for the study of medium 
modifications of nucleon properties. The nontopological soliton bag model, proposed by 
Friedberg and Lee, exhibits a dynamical bag formation due to the coupling of quarks to 
the phenomenological scalar field <fi. The quarks inside the soliton bag couple not only to 
the scalar field <fi that binds the quarks together into nucleons, but also to additional meson 
fields a, u, and p generated by the nuclear environment. The nucleons interact through the 
self-consistent exchange of a, u, and p meson fields that are treated as classical fields in the 
mean-field approximation. This model enables us to investigate the medium modification 
of nucleon properties because the soliton bag is significantly influenced by the additional 
meson fields in nuclear medium. 

We have considered two approaches for the cm. correction to nucleon properties so 
that we can estimate how sensitive the results are to the cm. correction approach used. 
We have adopted three parameter sets in the Friedberg-Lee model that are constrained 
by reproducing free nucleon properties. The quark-meson coupling constants are fitted to 
reproduce the empirical saturation properties of nuclear matter. The present model can 
provide a reasonable description of nuclear matter. We have found that the properties of 
the nucleon are significantly modified in nuclear medium. At normal nuclear matter density, 
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the nucleoli radius increases by about 10-16%, while the proton magnetic moment increases 
by about 11-16%. We have applied the present model to study the properties of spherical 
nuclei and found that it could give a reasonable description of the ground state properties 
of finite nuclei. It is very interesting to see the nucleon properties change from the nuclear 
surface to the nuclear interior. The present model incorporates explicit quark degrees of 
freedom into nuclear many-body systems. It is notable that the quark structure of the 
nucleon plays a crucial role in the description of nuclear matter and finite nuclei. 
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TABLE I: The nuclear matter properties in the present model with the parameter sets A, B, and 
C. The saturation density and the energy per particle are denoted by po and E/A, the symmetry 
energy by a sym , the incompressibility by K, and the effective mass by M*. 



Po E/A a sym K M*/M 
(fm" 3 ) (MeV) (MeV) (MeV) 



Set A 


0.15 


-16.0 


35 


302 


0.81 


Set B 


0.15 


-16.0 


35 


384 


0.72 


Set C 


0.15 


-16.0 


35 


184 


0.87 



TABLE II: The binding energy per nucleon E/A and the rms charge radius R c for 40 Ca, 90 Zr, and 
208 Pb. 





E/A (MeV) 


Rc (fm) 




Set A Set B Set C Expt. 


Set A Set B Set C Expt. 


40 Ca 


8.53 7.66 9.35 8.55 


3.42 3.46 3.38 3.45 


90 Zr 


8.36 7.81 8.97 8.71 


4.26 4.27 4.24 4.26 


208p b 


7.48 7.13 7.95 7.87 


5.52 5.50 5.54 5.50 



TABLE III: The single-particle energies of proton (p) and neutron (n) for 40 Ca. The experimental 
data are taken from Ref. [3^. All energies are in MeV. 





Set A 


Set B 


Set C 


Expt. 


Shell 


p n 


p n 


p n 


p n 


lsi/2 


35.6 43.7 


37.8 45.9 


34.7 42.8 


50±11 50.0 


l?3/2 


24.6 32.4 


25.7 33.5 


24.4 32.2 


34±6 30.0 


lPl/2 


23.5 31.3 


23.6 31.4 


23.7 31.6 


34±6 27.0 


l<4/2 


12.5 20.0 


12.8 20.3 


12.6 20.2 


15.5 21.9 


14/2 


10.4 18.0 


9.3 16.8 


11.4 19.0 


8.3 15.6 


2s l/2 


8.1 15.6 


7.9 15.2 


8.4 15.9 


10.9 18.2 
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TABLE IV: The spin-orbit splittings of proton (AE p ) and neutron (AE n ) for 40 Ca and 208 Pb. The 
experimental data are taken from Ref. [3a]. All quantities are in MeV. 

^Ca ^Pb 



AEp AE n AE„ AE n 

(ld 5/2 - ld 3/2 ) (14/2 ~ jj3/2) (lgg/g - 1^7/2) ( 2 /7/2 - 2/5/2) 
Set A 2.0 2.0 1.2 0.7 
Set B 3.5 3.5 2.2 1.3 
Set C 1.1 1.1 0.6 0.4 
Expt. 7.2 6.3 4.0 1.8 
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FIG. 1: The self- interaction of the soliton field as a function of 4> f° r the two parameter sets 
characterized by a = and B = 0. 
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FIG. 2: The cm. energy of the nucleon in nuclear matter, -E c .m. = E — M* , as a function of 
the effective quark mass, m* = g%o. The solid lines show the results of the first cm. correction 
approach where M* is given by Eq. (I30I) . The dashed lines show the results of the second cm. 
correction approach with M* given by Eq. (I32p . We note that m* = in free space and m* ~ —240 
MeV at normal nuclear matter density. 
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r(fm) 



FIG. 3: The quark wave functions (u, v) and the scalar field 4> as a function of the nucleon radius 
(r) with the three parameter sets. The results of the nucleon in free space are denoted by uq, vo 
and 4>q. The results at the density po and 2po are denoted by u\, v\, (pi and 112, V2, (p2, respectively, 
where po = 0.15 fm -3 is the normal nuclear matter density. 



1.0 
0.9 
^ 0.8 
0.7 
0.6 
0.5 

0.0 0.1 0.2 0.3 0.4 

p (fm- 3 ) 

FIG. 4: The ratio of the effective nucleon mass in nuclear matter to that in free space, M*/M, as 
a function of nuclear matter density p. The results with the parameter sets A, B, and C are shown 
by the solid, dashed, and dotted lines, respectively. 
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0.0 0.1 0.2 0.3 0.4 

p(fm- 3 ) 

FIG. 5: The ratio of the nucleon rms radius in nuclear matter to that in free space, r*/r c , as a 
function of nuclear matter density p. The lines are labeled as in Fig. HI 




p(fm- 3 ) 

FIG. 6: The ratio of the proton magnetic moment in nuclear matter to that in free space, p>*/p p , 
as a function of nuclear matter density p. The lines are labeled as in Fig. HJ 
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0.0 0.1 0.2 0.3 



p(fm- 3 ) 

FIG. 7: The ratio of the axial-vector to vector coupling constants in nuclear matter (<7a/<?v)* to 
that in free space gA/dv as a function of nuclear matter density p. The lines are labeled as in 
Fig.H 




FIG. 8: The charge density distributions for 208 Pb compared with the experimental data (solid 
line) [37]. The results with the parameter sets A, B, and C are shown by the dashed, dotted, and 
dot-dashed lines, respectively. 
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FIG. 9: The ratio of the proton rms radius in 208 Pb to that in free space as a function of radial 
coordinate r. The results with the parameter sets A, B, and C are shown by the solid, dashed, and 
dotted lines, respectively. 
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FIG. 10: The ratio of the proton magnetic moment in 208 Pb to that in free space as a function of 
radial coordinate r. The lines are labeled as in Fig. [9l 
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